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QCD with two flavours of massless colour-sextet quarks is considered as a model 

Oh| for conformal/walking Technicolor. If this theory possess an infrared fixed point, as 

(— ! ' 

indicated by 2-loop perturbation theory, it is a conformal(unparticle) field theory. 

C*~) ■ If, on the other hand, a chiral condensate forms on the weak-coupling side of this 

> , 

would-be fixed point, the theory remains confining. The only difference between 

00 
ON 



such a theory and regular QCD is that there is a range of momentum scales over 



CN ' which the coupling constant runs very slowly (walks). In this first analysis, we 

o ; 

simulate the lattice version of QCD with two flavours of staggered quarks at finite 
temperatures on lattices of temporal extent N = 4 and 6. The deconfinement and 

•l-H , 

^ ' chiral-symmetry restoration couplings give us a measure of the scales associated with 

confinement and chiral-symmetry breaking. We find that, in contrast to what is seen 
with fundamental quarks, these transition couplings are very different. f3 = Q/g 2 for 
each of these transitions increases significantly from Nt = 4 and Nt = 6 as expected 
for the finite temperature transitions of an asymptotically-free theory. This suggests 
a walking rather than a conformal behaviour, in contrast to what is observed with 
Wilson quarks. In contrast to what is found for fundamental quarks, the deconfined 
phase exhibits states in which the Polyakov loop is oriented in the directions of all 
three cube roots of unity. At very weak coupling the states with complex Polyakov 
loops undergo a transition to a state with a real, negative Polyakov loop. 
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I. INTRODUCTION 

The advent of the LHC era has led to renewed interest in extensions of the Standard Model 
which have a strongly-coupled symmetry-breaking (Higgs) sector. The most promising of 
these theories are the so-called Technicolor theories [jj, |2|. These are QCD-like theories 
with massless fermions (techni-quarks) in which the Goldstone bosons of spontaneously- 
broken chiral symmetry (techi-pions) play the role of the Higgs field, giving masses to the 
W and Z. Such models need to be extended in order to also give masses to the quarks and 
leptons. Phenomenological difficulties with such models, such as flavour- changing neutral 
currents which are too large, can be avoided if the fermion content of a candidate (extended) 
Technicolor theory is such that the running gauge coupling constant evolves very slowly. 
Such theories are referred to as Walking Technicolor models 3-6]. 

Let us consider a Yang-Mills gauge theory with Nf fermions in a specified (not-too-large) 
representation of the "colour" group. The evolution of the coupling constant g in such a 
theory is described by the Callan-Symanzik beta function /3(g) defined by 

dg g^ g^ 

^ )= ^ = - /3 °(4^- /3 W' W 

where /i is the momentum scale at which the running coupling constant g(fi) is defined. 
Pofli,... are given by perturbation theory. For Nf sufficiently small, (3q and (3\ are both pos- 
itive, the theory is asymptotically free and confining, and chiral symmetry is spontaneously 
broken. There exists some value of Nf above which /3 (and /3\) are negative and asymptotic 
freedom is lost. Between these two regimes is a range of Nf over which (3q is positive and 
f3\ is negative. In this range the theory remains asymptotically free, but if this 2-loop beta 
function describes the physics, it has a second zero which is an infrared (IR) fixed point and 
the theory is a conformal (unparticle) field theory. There is, however, another possibility. If 
the coupling becomes strong enough that a chiral condensate forms before this would-be IR 
fixed point can be reached, the fermions will become less effective at screening technicolor. 
f3 will then start to decrease again, and the theory will be confining in the infrared. There 
will, however, be a range of /i over which (3 is small and g evolves only slowly, i.e. it walks. 
Since the formation of a chiral condensate which spontaneously breaks chiral symmetry is a 
non-perturbative process, the boundary between conformal and walking behaviour cannot 
be determined perturbatively. Lattice gauge theory simulations can enable one to decide 
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between these two options for a theory with a specified gauge group, fermion technicolor 
representation and Nf. 

For SU(N) gauge theories, the most promising candidates for walking behaviour have 
fermions in the fundamental, adjoint, symmetric 2-index tensor or antisymmetric 2-index 
tensor representations of the gauge group. A good summary of what is known is given in 



reference 



7| • Estimates of the value of Nf below which a gauge theory walks and above which 



it is conformal have been made using various methods , none of which can be guaranteed 



to capture the full non-perturbative behaviour of QCD-like theories |8Hl3| . This has led 
people to use lattice gauge theory simulations to study this boundary. There have been a 
number of simulations of QCD with Nf fundamental quarks, with Nf large enough that 
conformal or walking behaviour might be expected 14H29j]. While progress has been made, 
the boundary value for Nf is still uncertain. Some simulations have been made of SU(2) 



gauge theory with 2 Dirac flavours of adjoint fermions 3CH37I]. While early indications are 



that this is a conformal field theory, there is still a great deal of uncertainty. 

For QCD with colour-sextet (symmetric tensor) quarks, (3i changes sign at Nf = 1^, 
and asymptotic freedom is lost at Nf = 3^j. Hence only Nf = 2 and Nf = 3 lie in the 
domain of interest. Nf = 3 is just below the value where asymptotic freedom is lost and is 
thus expected to be conformal. This leaves Nf = 2 as the most likely candidate for walking 
behaviour. DeGrand, Shamir and Svetitsky have simulated lattice QCD with 2 flavours of 
colour-sextet Wilson quarks 38M40| . Their initial results suggest that this is a conformal 



field theory. Fodor et al. have performed some preliminary simulations of lattice QCD with 
2 sextet quarks using domain- wall quarks [4 If . 

We are performing simulations of lattice QCD with 2 colour-sextet staggered quarks. 
Staggered quarks have the advantage over Wilson quarks in having a simple chiral order 
parameter, and no parameter tuning is needed to find the chiral limit. We are performing 
simulations at finite temperatures, where the deconfinement and chiral-symmetry restoration 
temperatures can be used as measures of the scales of confinement and chiral-symmetry 
breaking respectively. By varying the number of time slices N t we can determine whether 
these transitions scale as finite temperature transitions or whether they are bulk transitions. 
Preliminary results of these simulations were presented at Lattice 2009 {42 1. 

Our simulations indicate that, whereas for fundamental quarks the deconfinement and 
chiral-symmetry restoration transitions appear coincident, for colour-sextet quarks these 
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transitions are well separated. Chiral-symmetry restoration occurs at a much weaker cou- 
pling than deconfmement. This differs from what is seen with Wilson quarks by DeGrand, 
Shamir and Svetitsky 39] where the two transitions appear coincident. Such a separation 
has been reported in earlier simulations with colour-adjoint quarks 43|, l44| . and in SU(2) 
gauge theory with colour-adjoint quarks 45|]. Both transitions move to significantly weaker 
couplings when N t is increased from 4 to 6, which is what would be expected for finite 
temperature transitions governed by asymptotic freedom. This in turn favours the walking 
scenario. 

In the deconfined region, just above the deconfmement transition, we find 3 states, where 
the Wilson Line (Polyakov Loop) is oriented in the directions of the 3 cube roots of unity, 
similar to what occurs for quenched QCD or QCD with adjoint quarks. For N t = 4 only the 
state with a real positive Polyakov Loop appears stable. The other two states, while long- 
lived, appear to be only metastable. For N t = 6, all 3 states appear to be stable. Between 
the deconfmement and chiral-symmetry restoration transitions there is another transition 
where the 2 states with complex Polyakov Loops disorder into a state with a real, negative 
Polyakov Loop. Machtey and Svetitsky have argued that such additional states where the 
Polyakov Loop has arguments ±27r/3 and tt are to be expected, and have pres ented evidence 
for their existence and metastability in simulations using Wilson quarks [461 1 . 

In section 2 we discuss our simulations. Our results are described in section 3. Sec- 
tion 4 presents discussions and conclusions and identifies directions for ongoing and future 
investigation. 



II. LATTICE SIMULATIONS WITH SEXTET QUARKS 



For our simulations we use the simple Wilson gauge action 



1 - -Re(TrUUUU) 



(2) 



The fermion action is based on the unimproved staggered-quark action written formally as 
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where p=T,^ t]^ d ^ witn 

D^{x) = \[Uf\x)^{x + fi) - U™{x - Mix 



(3) 
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To allow tuning the number of flavours to values of Nj which are not multiples of 4, and to 
use a positive-definite operator for the transition to pseudofermions, this is replaced with 

Sf=ExH[P+m][-p+m]} N ^ X - (5) 

sites 

We use the RHMC algorithm for our simulations in which the fractional powers of the 
positive-definite Dirac operator are approximated to any desired accuracy by a rational 
approximation, and each trajectory is subjected to a global Metropolis accept/reject step, 
thus removing all dependence on the updating increment. 

We perform simulations on 8 3 x 4, 12 3 x 4 and 12 3 x 6 lattices, over a range of values of 
(3 = 6/g 2 which covers all transitions. For each lattice size we perform runs at quark mass 
m = 0.005, m = 0.01 and m = 0.02 in lattice units, to allow us to access the chiral limit. 
Away from the transitions our run lengths are typically 10,000 length- 1 trajectories for a 
given m and f3. Near deconfinement transitions run lengths of 50,000 to 200,000 trajectories 
are used for each (5 and m. Some runs of 50,000 trajectories have also been used near the 
transitions from negative Polyakov Loop states to complex Polyakov Loop states. 

We create deconfined states with positive Polyakov Loops by starting a run at /3 — 7.0 
(weak coupling) from a completely ordered state (all Us equal to the identity matrix). The 
configurations from these runs are used to start runs at progressively smaller /3s (and masses). 
The states with negative Polyakov Loops are obtained by starting a run at ft — 7.0 with all 
U s equal to the identity matrix, except for the timelike (7s on a single time slice, which are 
set to the matrix diag(l, —1, —1). 

The triplet Wilson Line (Polyakov Loop) is used to identify the position of the decon- 
finement transition. The chiral-symmetry restoring phase transition occurs at that j3 above 
which the chiral condensate vanishes in the chiral (m — > 0) limit. Since the chiral 

extrapolation is difficult to perform with the masses we use, we estimate the position of the 
chiral transition from the positions of the peaks in the chiral susceptibilities as functions 
of mass. 

X* = V [{(W) 2 ) - <W 2 ] (6) 

where the ipips in this formula are lattice averaged quantities and V is the space-time volume 
of the lattice. Since we use stochastic estimators for ipip, we obtain an unbiased estimator for 
this quantity by using several independent estimates for each configuration (5, in fact). Our 



6 



estimate of {ipip) 2 is then given by the average of the (10) estimates which are 'off diagonal' 
in the noise. 

III. RESULTS FROM SIMULATIONS 

A. N t = 4 

We perform simulations at a selection of values in the range 5.0 < < 7.0 on 8 3 x 4 
and 12 3 x 4 lattices. For each of the 3 chosen masses (0.005, 0.01, 0.02) the values of the 
Wilson Line and chiral condensate on the 8 3 x 4 and 12 3 x 4 lattices are so close that we 
conclude that finite size effects are negligible. Figure [1] shows the Wilson Line and chiral 
condensates as functions of (3 for each of the 3 masses (0.005, 0.01, 0.02) on a 12 3 x 4 lattice. 
In the deconfined phase, we have included only those states where the Wilson Line is real 
and positive, noting that runs which start in a state with a positive Wilson Line remain 
in this state for the duration of the run. We have not included the results for the 8 3 x 4 
lattice, since at the resolution of figured] the points for the 2 lattice sizes would be virtually 
indistinguishable. 

It is clear that the Wilson Line is very small below (3 pa 5.4, and rises rapidly shortly 
after this value for all 3 quark masses. This is taken as a signal for the deconfinement 
transition. It is also clear that chiral symmetry remains broken well beyond this point. 
Thus the deconfinement and chiral-symmetry restoration transitions are far apart, unlike 
what is observed for Wilson quarks, where they appear to be coincident [391 ] . 

Figure |2] shows a histogram of the Wilson Line at (3 = 5.42, m = 0.02, showing a clear 
two-state signal, suggesting that this (3 is very close to the transition. The separation of the 
two states suggests that this transition is a first-order phase transition. We conclude that 
at m = 0.02 the deconfinement transition is at (3 = (3d = 5.420(5). For m = 0.01, two-state 
signals are seen at (3 = 5.411 and (3 = 5.412 leading to an estimate (3d = 5.4115(5). Finally 
we note that for m = 0.005, (3 = 5.4 appears to lie below the transition while f3 = 5.41 
appears to be above the transition leading to an estimate (3d = 5.405(5). Thus the mass 
dependence of the deconfinement (3, (3d, is very weak. 

Now let us turn to the chiral transition. Because this is only expected to be a phase 
transition at m = 0, and the crossover becomes smoother as the quark mass is increased, 



7 



it is difficult to determine its position directly from the chiral condensate at the masses we 
use. This is made more difficult since it is clear from the measured condensates as functions 
of mass, that the mass dependence is far from linear, making it extrapolating to m = 
difficult if not impossible. We therefore use the peak in the chiral susceptibility defined in 
equation [6] as an estimate of the position of the crossover for finite mass. This is shown 
in figure [3] for the two smallest masses. From these graphs we estimate that the transition 
occurs at = 6.30(5) for both /3s. We thus estimate that the phase transition at m = 
occurs at j3 = j3 x — 6.3(1). Note that the spacing of the /3s in this range is too large to allow 
us to use Ferrenberg-Swendsen reweighting to determine the positions of these transitions 
with more resolution. (The distributions of values of the plaquette action for adjacent (3s 
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FIG. 1: Wilson line (Polyakov Loop) and (ipip) as functions of /3 on a 12 3 x 4 lattice. 
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FIG. 2: Histogram of the Wilson Line at f3 = 5.42, m = 0.02 on a 12 3 x 4 lattice. 

do not overlap in this region.) 

At each of our two larger masses, we have performed a series of runs starting from 
(3 = 7.0 with a negative Wilson Line. From (3 = 7.0 down to (3 = 6.0, the Wilson Lines 
remain negative over the length of the 10,000 trajectory run for each ((3,m). By (3 = 5.8, 
these states have transitioned to states in which the Wilson Line is oriented in the direction 
of one of the complex cube roots of unity. Hence we deduce that there is a transition at 
(3 w 5.9. More discussion of this transition is to be found in the N t = 6 subsection. On 
the 12 3 x 4 lattice with m = 0.02 these states with complex Wilson Lines persist down to 
(3 = 5.48, and appear stable over at least 50,000 trajectories. For m = 0.01 these persist 
down to /3 = 5.46. For f3 < 5.46 at m = 0.02 or f3 < 5.45 at m = 0.01 but above the 
deconfinement transition, these states with complex Wilson Lines appear to be metastable 
and eventually decay into states with positive Wilson Lines. Figure |4] shows an example of 
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FIG. 3: Chiral susceptibilities as functions of /3 on a 12 3 x 4 lattice for m = 0.005,0.01, for a 
/3 range which includes the chiral transition. 



such metastability in our 12 3 x 4 simulations. As is to be expected, this metastability starts 
at larger values on an 8 3 x 4 lattice. We have observed no cases where configurations with 
positive Wilson Lines evolve to configurations with complex Wilson Lines for /3 values above 
the deconfinement transition. 



B. N t = 6 

We perform simulations on a 12 3 x 6 lattice at quark masses m = 0.005, m = 0.01 and 
m = 0.02 for values of /3 = Q/g 2 in the range 5.0 < < 7.2. We perform two series of 
runs starting at = 7.0, for m = 0.01,0.02. The first starts with the Wilson Line real 
and positive, and the second with the Wilson Line negative. For the lowest quark mass 
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FIG. 4: Time evolution of the argument of the Wilson Line (Polyakov Loop) at m = 0.02, /3 = 5.45 
showing a decay of a state with a complex Wilson Line (argument close to 2n/3) to a state with a 
real positive Wilson Line (argument close to 0). 

m = 0.005 we only perform one set of runs starting with a positive Wilson Line at /3 = 7.0. 

Above the deconfinement transition we again see a 3-state signal where the Wilson Lines 
orient themselves in the directions of one of the 3 cube roots of unity. A scatterplot showing 
such a 3-state signal is shown in figure [5^. Unlike the N t = 4 case, there is no sign of 
metastability and transitions between all 3 states are seen over the duration of our runs, up 
to (3 values far enough above the deconfinement transition that the mean relaxation time 
between tunnelings exceeds the lengths of our runs (50,000 trajectories). Figure [5b shows an 
example of such tunnelings. Thus to make meaningful measurements of the Wilson Line, we 
need to separate these 3 states. This we do by binning the Wilson Lines measured at the end 
of each trajectory according to their phase <fi into bins — 7r < < — n/3, —tt/3 < <fi < 7r/3, 
7r/3<</><7r. To increase our statistics, we use symmetry to include the complex conjugates 
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FIG. 5: a) Scatterplot of Wilson Lines for m = 0.02 and j3 = 5.58 on a 12 3 x 6 lattice showing the 
3-state signal. 

b) 'Time' evolution of the argument (phase) of the Wilson Lines for one of the 2 runs included in 
part (a). 
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of the Wilson Lines in the first of these bins, with the Wilson Lines in the last of these bins. 
Other observables are binned according to the values of the corresponding Wilson Lines. 

Figure Eh shows the Wilson Lines (Polyakov Loops) and chiral condensates as functions 
of (3 for all 3 masses from our simulations on a 12 3 x 6 lattice for the states with positive 
Wilson Lines (— 7r/3 < <f) < tt/3). Figure [6b shows the magnitudes of the average Wilson 
Lines and the chiral condensates for the states with complex or negative Wilson Lines. 
The deconfinement transition manifests itself as a rapid increase in the Wilson Line, which 
is clearly seen at a just above 5.5. The fact that the Z 3 centre symmetry is broken 
manifests itself in the difference between the magnitudes of the Wilson Lines in the positive 
Wilson Line and complex Wilson Line states. The rapid change in the magnitude of the 
complex/negative Wilson line between (3 = 6.2 and (3 = 6.6 marks the transition between 
a state whose Wilson Line is oriented in the direction of one of the complex cube roots of 
unity and one where the Wilson Line is real and negative. The chiral transition, above which 
(ipip) vanishes in the chiral (m — > 0) limit, clearly occurs at a /3 appreciably larger (3d- 

To determine the position of the deconfinement transition we examine histograms of the 
magnitudes of the Wilson Line close to the transition for each mass. Although the explicit 
breaking of the Z 3 centre symmetry means that the magnitudes of the positive and complex 
Wilson Lines are not identical, they are close enough in the vicinity of the deconfinement 
transition that this fact can be ignored. Such histograms are presented in figure [7J For 
the lower 2 masses the histogram for each beta represents 50,000 trajectories. At m — 0.02 
the histograms for (3 = 5.55 and (3 = 5.57 represent 100,000 trajectories each, while that at 
(3 = 5.56 is for 200,000 trajectories. In each case the histogram peaks at a low value below the 
transition and an appreciably higher value above the transition. Very close to the transition 
the peak of the histogram is relatively flat, with some hint of a double peaked structure. 
From these graphs we estimate that the transition (3s {(3d) are (3d(m = 0.005) = 5.545(5), 
(3 d (m = 0.01) = 5.550(5) and (3 d (m = 0.02) = 5.560(5), respectively. As in the N t = 4 case, 
the mass dependence is weak. 

Again we estimate the position of the chiral-symmetry restoration transition by examining 
the peaks in the chiral susceptabilities for m = 0.005, 0.01. As for N t = 4 we obtain estimates 
of the chiral condensate from 5 independent noise vectors, which yields an unbiased estimate 
for the chiral susceptibility. These chiral susceptibilities are plotted versus (3 in figure [8] for 
each of these lowest two masses. Again, the /3s we use in this region are two sparse to allow 
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FIG. 6: a) Wilson Line and chiral condensate for the state with a real positive Wilson Line as 
functions of f3 for each of the 3 masses on a 12 3 x 6 lattice. 

b) Magnitude of the Wilson Line and chiral condensate for the states with a complex or a real 
negative Wilson Line as functions of for each of the 3 masses on a 12 3 x 6 lattice. 



the use of Ferrenberg-Swendsen reweighting to interpolate between them. Since our estimate 
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FIG. 7: Histograms of the magnitudes of Wilson Lines for j3 values bracketing the deconfinement 
transition on a 12 3 x 6 lattice for a) m = 0.005, b) m = 0.01, c) m = 0.02. 

for the peak of the susceptibility plots for each mass is /3 = 6.60(5), we estimate for the 
position of the chiral transition at m = is /3 = (3 X = 6.6(1). Note that this estimate is for 
the states with positive real Wilson Lines only. 

Starting at (5 = 7.0, m = 0.02 with a real negative Wilson Line, we find that this state 
is stable for runs of up to 50,000 trajectories for /3s down to /3 — 6.6. Decreasing this to 
(5 = 6.4, we find that the system has transitioned to a state with a Wilson Line oriented in 
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FIG. 8: Chiral susceptibilities Xm> as functions of (3 on a 12 3 x 6 lattice for m = 0.005,0.01, for a 
/3 range which includes the chiral transition. 

the direction of one of the complex cube roots of unity. Figure [9] shows the time evolution 
of the phases of the Wilson Lines at = 6.4 and /3 = 6.6 for m = 0.02. We see that the way 
the transition proceeds is that the fluctuations in the phase become large as /3 approaches 
the transition from above until eventually they reach 27r/3 and 4n/3(— 2n/3). When this 
occurs the system spends appreciably more time at these two values that at intermediate 
values. We can then consider the system as being in a state with its Wilson Line oriented in 
the direction of one of the two complex cube roots of unity. On the finite lattice it tunnels 
between these 2 states. Approaching this transition from below we could describe it as the 
disordering of the 2 states with Wilson Lines in the directions of the complex cube roots of 
unity, as suggested by Machtey and Svetitsky. The position of this transition for m = 0.02 
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FIG. 9: The 'time' evolution of the argument (phase) of the Wilson Line for m = 0.02 on a 12 3 x 6 
lattice, a) for f3 = 6.4 and b) for f3 = 6.6. The horizontal lines are at 27r/3, tt and 4ir/3. 

is at (3 ~ 6.5, while for m = 0.01 it occurs at /3 ~ 6.4. 
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IV. DISCUSSION AND CONCLUSIONS 

We simulated QCD with two flavours of staggered quarks at finite temperatures using the 
RHMC method. Our simulations were performed on 8 3 x 4, 12 3 x 4 and 12 3 x 6 lattices with 
3 quark masses, to allow for chiral extrapolations. We find widely separated deconfinement 
and chiral-symmetry restoration transitions. Both the deconfinement and chiral transitions 
move to significantly lower couplings as N t is increased from 4 to 6, which is the expected 
behaviour for finite temperature transitions in an asymptotically free theory. This suggests 
that the theory is confining with spontaneously-broken chiral symmetry, while being under 
the control of the weak-coupling asymptotically-free ultraviolet fixed point, i.e. that it walks. 

The deconfinement transition occurs at /3 = fl d where for Nt = 4, f3 d {m = 0.02) = 
5.420(5), p d (m = 0.01) = 5.4115(5) and (3 d (m = 0.005) = 5.405(5) while for N t = 6, 
p d (m = 0.02) = 5.560(5), d {m = 0.01) = 5.550(5) and (3 d (m = 0.005) = 5.545(5). The 
chiral-symmetry restoration transition occurs at = (3 X , where /3 X = 6.3(1) at N t = 4 and 
X = 6.6(1) at N t = 6. 

The large separation of the deconfinement and chiral transitions indicates that the en- 
hanced attraction between quark- ant iquark pairs over that for fundamental quarks due to 
the larger quadratic Casimir operator for the sextet representation (10/3 versus 4/3), causes 
a chiral condensate to form at a distance much shorter than the scale of confinement. This 
effectively removes the quarks from consideration at longer distances where the theory will 
behave more like a pure glue (quenched) theory. This is presumably why, in the deconfined 
phase, we see a three state system, where the Wilson Lines align in the directions of one of 
the cube roots of unity, a relic of the now-broken Z3 symmetry. The breaking of Z3 is seen 
in the difference in magnitudes of the real positive Wilson Lines versus those with phases 
close to ±27r/3. At N t = 4, the states with complex Wilson Lines are only metastable, while 
at Nt — 6 all 3 states appear to be stable. The existence of states with all 3 Z s Wilson Line 



orientations has been predicted by Machtey and Svetitsky 46] who observed them in their 
simulations with 2 flavours of colour-sextet Wilson quarks. They also observed the metasta- 
bility of the states with complex Wilson Lines. Earlier simulations by DeGrand, Shamir and 
Svetitsky had reported deconfined states with Wilson Lines oriented in the directions of the 
complex cube roots of unity 



More work needs to be done to determine whether the chiral-symmetry restoring transi- 
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tion of these complex/negative Wilson Line states is coincident with that of the state with 
a positive Wilson Line. In addition we would like to know whether this chiral transition 
shows the scaling properties of the expected 0(2)/ 0(4) universality class. The fact that this 
transition occurs at relatively weak coupling should help us in this endeavour. 

We have also observed an additional transition. At /3 ~ 5.9 on an N t = 4 lattice with 
m = 0.02 or m = 0.01, the states with complex Wilson Lines disorder to a state with a 
negative Wilson Line. Such a transition is also observed for N t — 6 at /3 ~ 6.5 for m = 0.02 
and ~ 6.4 for m = 0.01. The existence of states with negative Wilson Lines and of such 
a transition is predicted and observed by Machtey and Svetitsky. Such a transition would 
be expected to be either a second-order phase transition in the universality class of the 
3-dimensional Ising model, or a first-order phase transition. The large increase in the /3 for 
this transition between N t = 4 and N t — 6 makes us suspect that this is a lattice artifact. 
We also note that by N t = 6, it is close to the chiral transition and could well merge with 
it at larger ]V t . Since the negative Wilson Line state (phase it) comes from disordering 
the two states with phases ±27r/3 (a fact also predicted by Machtey and Svetitsky), the 
magnitude of the Wilson Line above the transition is approximately half what it is below 
the transition. Just below the transition, the magnitude of Wilson Line in the complex 
Wilson Line states is approximately 2/3 that of the positive Wilson Line state, so that after 
the transition the magnitude of the negative Wilson Line is approximately 1/3 of that of 
the positive Wilson Line. This would suggest that the transition might be associated with 
the symmetry breaking SU(3) ->• SU(2) x U(l). 

We need to understand why we see well-separated deconfinement and chiral transitions 
with staggered fermions, whereas DeGrand, Shamir and Svetitsky observed these transitions 
to be coincident with Wilson fermions. Of course, since we are far from the weak-coupling 
limit, different fermion actions do not have to have the same physics. This is especially 
true, if we happen to be in the strong-coupling domain, beyond the infrared fixed point of 
a conformal field theory. Of course, the observation from our simulations that the coupling 
at each of the transitions is decreasing as the lattice spacing is decreasing would appear 
to exclude this possibility, since the /3 function becomes positive above this fixed point. 
However, it has recently been suggested that there could be two non-trivial fixed points {^tJ 
in such theories. In this case, if we are beyond the second non-trivial fixed point (which would 
be an ultraviolet fixed point), the coupling would decrease at short distances. Of course, 
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if we are beyond the region where the ultraviolet behaviour of the theory is controlled by 
asymptotic freedom, drawing any conclusions is pure speculation. 

To better understand our results, and to help determine whether the theory is indeed 
walking, rather than conformal as indicated by the work of DeGrand, Shamir and Svetitsky 
using Wilson fermions, we are now extending our simulations to N t = 8, where finite lattice 
spacing effects should be reduced, to see if the deconfinement and chiral-symmetry restora- 
tion transitions remain consistent with being finite-temperature transitions. We also plan to 
perform simulations at zero temperature, measuring the chiral condensate, the hadron spec- 
trum, f n , etc., to test other aspects of the theory which should help us determine whether 
this theory is conformal or walking. In addition we will try to determine the running of a 
suitably-defined renormalized coupling. Having two different spatial lattice sizes at Nt = 4 
showed us that finite size effects are small. We need a second spatial lattice size for N t = 6. 
We have learned that the authors of reference [4l| are now starting simulations of this theory 



using improved staggered fermions, which should help resolve these issues 48]. 

We are now performing simulations of lattice QCD with 3 flavours of staggered quarks 
at finite temperature. Since this theory is almost certainly conformal, it is interesting to 
determine whether its behaviour is qualitatively different from that of the Nf = 2 case, and 
whether our simulations can see this conformality. 
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